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ABSTRACT 

Interpretation of imagery of the solar chromosphere in the widely used Ca II 854.2 nm infrared 
line is hampered by its complex, three-dimensional and non-LTE formation. Forward modelling is 
required to aid understanding. We use a 3D non-LTE radiative transfer code to compute synthetic 
Ca II 854.2 nm images from a radiation-MHD simulation of the solar atmosphere spanning from the 
convection zone to the corona. We compare the simulation with observations obtained with the CRISP 
filter at the Swedish 1-m Solar Telescope. We find that the simulation reproduces dark patches in the 
blue line wing caused by Doppler shifts, brightenings in the line core caused by upward- propagating 
shocks and thin dark elongated structures in the line core that form the interface between upward 
and downward gas motion in the chromosphere. The synthetic line core is narrower than the observed 
one, indicating that the sun exhibits both more vigorous large-scale dynamics as well as small scale 
motions that are not resolved within the simulation, presumably owing to a lack of spatial resolution. 
Subject headings: methods: numerical — Sun: atmosphere — Sun: chromosphere 



1. INTRODUCTION 

The solar chromosphere has traditionally been ob- 
served in the Ca II H&K and the Ha lines. The use of 
these lines as a diagnostic suffers from significant draw- 
backs. The Ca II H&K lines have wavelengths in the 
violet part of the spectrum. A trade-off has to be made 
between filter width and exposure time due to the lack 
of photons and the lower sensitivity of CCD chips at 
these wavelengths. Either the filter has to be wide, lead- 
ing to significant low-chromospheric line wing contribu- 
tions, or the long exposure time leads to lower spatial 
resol ution due to difficulties with the im age restoration 
(e.g. iLites et~aT1 119991 : IWoger et all 12006( 1. In addition, 
the Ca II H&K lines are subject to partial redistribu- 
tion (PRD), which complicates proper modeling of these 
lines. 

Observations in the Ha line suffer to a much lower ex- 
tent from these effects, but Ha line formation is very 
complex and proper modeling requires non-equilibrium 
ionization of hydrogen to be taken into account, physics 
that so far has only be en treated in detail i n one- 
dimensional (ID) models ((Carlsson fc Steinll2002(l . 

These problems are much less severe for the 
Ca II 854.2 nm line. Non-equilibriu m and PRD effect s 
are less important for its formation (|Uitenbroekl Il989(l . 
Observations can be done with narrow filters and short 
exposure times, yielding a clean and high resolution view 
of the chromosphere. This makes it an excellent diagnos- 
tic, both from the observational as the modeling point-of- 
use, albeit at a lower diffraction-limited resolution than 
the Ca II H&K and the Ha lines. 

Fabry-Per ot instruments such a s IBIS (|Cavalrini| [2006) 
and CRISP (|Scharmer et aLlfeOOSD now routinely observe 
the Ca II 854.2 nm line with high spectral, spatial, and 
temporal resolution. However, interpretation of data 
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obtained with these instruments is non-trivial; forward 
modeling by computation of synthetic filtergrams from 
radiation-MHD simulation is needed to gain insight in 
the formation of the line. 

The for mation of th e Ca I I 854.2 nm line has been 
studied by lUitenbroekl (120061) in t he ID semi-empirical 
FALC mo del dFontenla et al.lll993l) and the ID dynamic 
models of lCarlsson fc Steinl (Il999f) a nd a 3D snapshot of 
convection by lAsplund et alj (|2000D . The 3D model did 
not include magnetic fields and extended up to 1 Mm 
above T500 = 1, lower than the formation height of the 
line core in the FALC model (see Fig. 5 of E 

auzzi et al.l 

12008(1 . In this paper we present synthetic images of the 
Ca II 854.2 nm line computed with a 3D non-LTE radia- 
tive transfer code from a snapshot of a radiation-MHD 
simulation that extends from the convection zone to the 
corona, capturing the whole height range of formation of 
the line. 

2. OBSERVATIONS 

Figure [3] shows observations in the Ca II 854.2 nm 
line obtained with the C RISP imaging spectropolarime- 
ter ((Scharmer et al.ll2008H at the Swedish 1-m Solar Tele- 
scope. The data was obtained on June 13, 2008 and was 
taken at disk-center in a coronal hole. A scan through 
the line was made using 24 filter positions from —0.92 
to +0.19 A with a FWHM of the filter of 110 mA. The 
image quality was imp roved by post-processin g using the 
MOMFBD algorithm ((van Noort et al.ll2005h . 

3. RADIATIVE TRANSFER COMPUTATION 

The radiative transfer computation was based on 
a snapshot compute d with the Oslo Stagger Code 
((Hansteen et al.ll2007( l . This code employs an LTE equa- 
tion of state and includes non-LTE radiative losses using 
a multi-group opacity method and thin radiative cool- 
ing in the corona and upper chromosphere. It includes 
thermal conduction along magnetic field lines. The sim- 
ulation had 256 x 128 x 160 grid points, corresponding 
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Fig. 1. — Average profile of the Ca 8542 IR line. Solid red: 
FTS atlas; solid black: simulation. The simulated continuum has 
a lower intensity and the line core is narrower than in the observed 
profile. 

to a physical size of 16.6 x 8.3 x 15.5 Mm. The snap- 
shot has an averaged unsigned magnetic field strength 
of 150 G and covers a height range from 1.5 Mm be- 
low the photosphere up to 14 Mm above it, spanning 
from the upper convection zone up into the corona. For 
the radiative transfer computation the grid was inter- 
polated onto a grid with 213 points in the z-direction 
covering a height range between -0.5 and 5.3 Mm rela- 
tive to T500 = 0. The Ca II model atoms consists of 5 
bound levels plus a continuum with 5 lines (Ca II H&K 
and the infrared triplet) and 5 bound-free transitions. 
The Ca II 854.2 nm line had 200 frequency points. All 
lines were computed in complete redistribution. Velocity 
fields were taken into account and no additional micro- 
turbulence was added. The electron density was com- 
puted assuming LTE ionization for all relevant species. 
Photoionization by hydrogen Lyman lines was not taken 
into account. The ray quadra ture was computed using 
the A2 set of ICarlsor] (|1963f h The radiative transfer 
computation was performed using an MPI-parallelized, 
domain-decomposed code named MULTI3D. This code 
is ba sed on the ID code MULTI (|Scharmer fc Carlssonl 
119851 ) and includes the same physics, but has a 3D short- 
characteristics formal solver, allowing evaluation of the 
3D radiation field. 

4. RESULTS 

Figure [T] compares the average quiet s un profile of the 
Call 854.2 nm line from the FTS atlas (|Neckel k. Labs! 
fl98l to the average profile computed from the simula- 
tion. The simulated continuum is 3.5% less bright than 
the observations, the line wings fit rather well, while the 
line core is narrower and the bisector does not show the 
inverse-C shape of the observations. 

The vertically emergent intensity in the Ca II 854.2 nm 
line was convolved with the CRISP filter transmission 
function. Figure [2] compares the observed and synthetic 
profile averaged over the simulation box and a quiet 
region of the same size in the observations (shown in 
Fig. [3]). The synthetic profile has a less steep inner wing 
and a narrower line core. 

Figure [3] compares a quiet sun region from the observa- 
tions to the synthetic filtergrams. The left panels show 
reversed granulation with several magnetic bright points. 
The contribution function in the upper-left panel of Fig. 0] 
shows that the intensity in the magnetic bright points 
are formed deeper than in the surrounding low-field- 




Fig. 2. — Line core of the Ca 8542 IR line. The profiles have been 
normalized on their values at AA = —0.92 A. Solid red: CRISP; 
solid black: simulation; plus signs indicate the different filter posi- 
tions. The gray curves show the shape and positions of the filter 
used to produce the synthetic images in Fig. [3] The synthetic line 
core is narrower than the observed one. 

strength atmosphere , similar to G band bright points 
(jCarlsson et al.ll2004h . 

The middle panels of Fig. [3] show the atmosphere as 
seen in the "knee" of the line, where the profile changes 
from the LTE win gs to the non-LTE formed line core 
(|Cauzzi et al.ll206"8[) . Because of the difference in the line 
widths in the simulations and the observations this is at 
different offsets relative to line center, at —0.39 A for the 
observations and —0.14 A for the simulations. Neverthe- 
less, they show a remarkably similar scene. The bright 
background shows the upper photosphere and lower chro- 
mosphere at around z = 0.5Mm where the reversed gran- 
ulation pattern ch anges to one dominated by acoustic 
shock interference (|Wedemever et aLll2004lh 

The upper-middle panel of Fig. [3] shows black struc- 
tures, for example at (15, 1) and at (15.5, 13) Mm. Sim- 
ilar roundish or thin elongated structures appear in the 
synthetic image. The corresponding contribution func- 
tion panel in Fig.|4]shows that they are formed at around 
z = 1.5 Mm, 1 Mm higher than the bright background. 
Comparison with the lower right panel shows that the 
dark patches are associated with upflows and the line 
profile in Fig. [5] shows a blueshifted line at the positions 
of the dark patches. 

There are three brightenings along the cut in the lower 
middle panel of Fig. [3j at x = l, x — 9.5 and a; = 15.5 Mm, 
all located above a photosphcric magnetic field concen- 
tration (middle right panel of Fig. [4j) . 

The right panels show the line core. The simulated 
image has a higher contrast, which can at least partly 
be explained by the fact that the observations are not 
corrected for scattered light. The observations show an 
amorphous background with a number of brightenings, 
such as at (x,y) = (3.5,3.5) and (16, 8) Mm. They are 
often, but not always, located on top of or near bright 
points in the photosphere. The simulated line core im- 
age (lower right) shows similar brightenings, for example 
at (x,y) — (1,3) Mm. Fig. [4] shows that this brighten- 
ing is associated with a strong downflow above a mag- 
netic element. This brightening is formed much deeper 
than is the intensity along the rest of the cut. A simi- 
lar brightening at (a;, y) = (9.5, 1.5) Mm is not associated 
with a photosphcric magnetic field concentration, how- 
ever. The corresponding line profile in Fig. [5] shows an 
emission peak slightly bluewards of line-center and an 
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Fig. 3. — Observed (top row) and synthetic (bottom row) Ca II 854.2 nm images at different positions in the line. Left: at AA = — 0.87 A; 
middle: in the "knee" of the line, for the observations at AA = — 0.39 A and for the synthetic image at AA = — 0.14 A; right: in the line 
core at AA = oA. The brightness scales are clipped at the brightest and darkest 0.1% of the pixels. The observed and synthetic images 
show similar structure in all panel pairs. The wavelengths of the middle panel pair are in the "knee" of the line profile, and because of the 
difference in the line core width this is not at the same wavelength. The red line indicates the cut displayed in Figs. [4] and [5] 
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Fig. 4. — Vertical cut through the atmosphere along the red line in Fig. [3] The left-hand panels show the contribution function to the 
intensity within the CRISP filter, on an inverted brightness scale and with each column scaled to maximum contrast for improved visibility. 
The overplotted grey curves indicate the emergent intensity in arbitrary units. The right-hand panels show, from top to bottom, the gas 
temperature, the absolute magnetic field strength and the vertical velocity (positive means upfiow). The red curve in the velocity panel 
indicates the zero vertical velocity contour. 
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Fig. 5. — Emergent spectrum of the Ca II 854.2 nm line along 
the cut indicated in Fig. [3] with positive velocity corresponding to 
a redshift. The brightening at (x, y) = (1, 3) Mm in the lower right 
panel of Fig.[3]is associated with a 12km/s redshift of the line core. 
The dark patches between y = 5 and y = 10 Mm along the cut in 
the lower middle panel of Fig. \3\ are associated with a blueshift of 
5 km/s 

intensity dip on the red side of the line typical of an up- 
ward p ropagating wave with d ownflowing material above 
it (see iCarlsson fc Steinlll997l for a detailed description 
in the case of the Ca II H2V grains). These shocks can be 
either acoustic or magneto-acoustic, the radiative trans- 
fer mechanism causing them to be bright in the line core 
remains the same, and a single snapshot alone is not suf- 
ficient to distinguish between them. Analysis of a time- 
series is needed to determine the wave type. 

On top of the background in the upper right panel 
of Fig. [3] one can see filamentary dark structures that 
can be both straight and curved. They are mainly con- 
centrated in the top left part of the image, but can 
be found throughout the field of view. The simulated 
line-core image shows similar dark filamentary struc- 
tures (for example at (x,y) = (2.5,3) and (10, 6) Mm). 
They appear similar to the observed ones, but tend to 
be straighter, though a slightly curved example can be 
found at (at, y) = (4, 5) Mm. They are the interface be- 
tween a patch of upflowing and downflowing material. 
They form where the vertical velocity at their formation 
height is zero and the line core is located at the rest fre- 
quency of the line. On either side of the dark filament 
the gas is moving up or down, shifting the line core to 
the blue or red and thus showing the brighter line wings 
at the rest frequency. 

5. DISCUSSION & CONCLUSIONS 

The simulation reproduces a number of observed fea- 
tures rather well, but there are a number important dif- 
ferences that carry large diagostic potential. 

The difference in the continuum intensities in Fig.[T]can 
be explained by a small difference in the effective temper- 
ature of the simulation and the sun combined with the 



effect of box oscillations in the simulation, which causes 
fluctuations in the continuum intensity with an ampli- 
tude of about 1%. 

The simulated spatially-averaged line-core profile is 
narrower than the observed one (Fig. [21). We computed 
the deviations of the position of the line core from the 
rest-wavelength over the field-of-view for both the simu- 
lation and the observation. The standard deviation for 
the simulations is 1.1 km/s, which is 50% smaller than for 
the observations. This indicates that the solar chromo- 
sphere has more vigorous large-scale dynamics than the 
simulation. However, this effect is too small to explain 
the difference in average line-core width. 

We also inspected the line-core widths for individual 
resolution elements. The simulated line cores are on av- 
erage narrower than the observed ones. This is the dom- 
inant reason for the difference in the average core width, 
indicating that the sun is dynamic on scales not resolved 
by the simulation. 

Close inspection of the line core images of Fig. [3] shows 
that the observations have more structure on scales of 
about 0.1 Mm than the simulation, also hinting that the 
simulation does not have sufficient resolution to generate 
enough small-scale motion to broaden the line core. 

In summary, we have computed the 3D, non-LTE, 
emergent intensity in the Ca II 854.2 nm line from a 
snapshot of a radiation-MHD simulation. It reproduces 
a number of observed features rather well, i.e., black 
patches in the blue 'knee' of the line, and thin dark lines 
and brightness enhancements in the line core. However, 
the simulated line core is narrower than the observed one, 
indicating that the sun shows more vigorous large scale 
dynamics and motions at smaller scales than can be re- 
solved in the simulation. We intend to repeat the same 
analysis with a radiation-MHD simulation with higher 
spatial resolution. 
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